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Abstract 

We obtain an analytic expression for the specific heat of a system of N rigid 
rotators exactly in the high temperature limit, and via a pertubative approach in 
the low temperature limit. We then evaluate the specific heat of a diatomic gas with 
both translational and rotational degrees of freedom, and conclude that there is a 
mixing between the translational and rotational degrees of freedom in nonextensive 
statistics. 
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I Introduction 



The nonextensive generalization of Boltzmann-Gibbs statistics was proposed by Tsallis 
[lj to study systems involving long-range interactions, long-range microscopic memory, 
and nonequilibrium phenomenon. The nonextensive statistics has been applied in a wide 
range of areas like self-gravitating systems [21 El H], solar neutrinos [SJ Ej, and biological 
systems pD, [8]. The functional form of the generalized entropy reads 

= kln q W, q £ R+, (1.1) 

where k is a constant generalizing the Boltzmann constant k B for an arbitrary q, and W 
is the weight factor. In the q — > 1 limit the Tsallis entropy S q reduces to the standard 
Boltzmann-Gibbs entropy S. It is interesting to observe that these two entropies may be 
interrelated via hypergeometric functions as 

S/k B = 2^(1, 1; 2; -(1 - q)SJk) S q /k, (1.2) 
Sjk = iF x {l\ 2; (1 - q)S/k B ) S/k B . (1.3) 

Based on the entropy (11.11) a canonical formulation was developed in PUD]- As a concrete 
example of using this formalism for studies of various statistical systems at equilibrium 
the classical monoatomic ideal gas has been investigated as a iV-body problem [TT1 [T2] in 
detail, and, in particular, the specific heat of this system has been computed exactly. 

In the above context study of the rigid rotator has been initiated in [T3], HI] . Using 
molecular dynamics a numerical study of iV rotators with a long-range interaction has 
been done in [13], where a long-standing metastable state with anomalous duration sug- 
gesting nonextensivity has been observed. In [H] an analytical study of the rigid and 
non-rigid rotators has been done. These authors, however, obtained the specific heat of 
the system by studying a single rigid rotator. In the regime of nonextensive statistics, 
strictly speaking, the study of a single rigid rotator cannot be extended to a system of 
N rigid rotators as it is familiarly done in the extensive Boltzmann-Gibbs scenario. A 
systematic study of statistics of N rigid rotators is, therefore, essential for understanding 
the dependence of the specific heat on the temperature. We will comment on the other 
relevant issues in appropriate places. 

In the present work we first evaluate the specific heat at the high temperature limit 
for N rigid rotators, and a diatomic gas endowed with both translational and rotational 
degrees of freedom. Towards this end we construct an implicit equation involving the 
generalized partition function and the internal energy. In the high temperature limit the 
implicit equation is solved exactly and explicit expression for the specific heat is evaluated. 
Contrarily, obtaining an exact solution of the said implicit equation for a system of N rigid 
rotators in the low temperature limit appears to be difficult. We, therefore, use the method 
developed in [10] that interrelates the ensemble probabilities obtained via the second and 
the third constraint approaches, respectively, after performing a suitable transformation 
on the temperature. This transformation connects the second constraint approach [9] that 
utilizes unnormalized g-expectation values to the third constraint [10] approach based on 
the normalized g-expectation values. To construct the pertinent transformation here we 
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adopt a pertubative approach by disentangling the g-exponential in an infinite series of 
ordinary exponential [15], and retain terms upto the order (1 — q) 2 . In the context of the 
gas of N diatomic molecules we notice that a mixing of the translational and rotational 
degrees of freedom in the nonextensive terms occur for the specific heat. We here assume 
a large mass for the diatomic molecules so that their vibratory motions may be neglected. 
The plan of the paper is as follows. The specific heat of N rigid rotators in the high and low 
temperature regimes are discussed in Sections [III and IIII[ respectively. This is followed by 
our discussions on the specific heat of a diatomic gas in Sec. [TV] Our concluding remarks 
are given in Sec. [V] 



II Rigid rotators in the high temperature limit 

The energy eigenvalues of an ensemble of N isotropic rigid rotators read 
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where I is the moment of inertia of a rotator. The generalized partition function of 
a canonical ensemble of N rotators at an inverse temperature (3 = 1/kT in the third 
constraint framework [10] is given by 



l-(l-q)^(C-n R U t 



1-q 



c = £[pf(/?)] 9 , (2-2) 



where Qr = 2I/h 2 ,j3 = (3/£Ir. The degeneracy factor T> and the sum C read 
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In the expression (12.21) the ensemble probability in the third constraint approach is given 
by 



i-(i-q)^(c-n R u q ) 
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At high temperature limit the rotational levels are closely spaced, and, therefore, the sum 
in (12.21) is well-approximated by an integral 
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The integral may be exactly evaluated as 
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where A( fc ) = 1 + k(l — q)N. As a consequence of the property [TU] 



(2.7) 



the above expression of the generalized partition function may be regarded as an implicit 
equation. The internal energy of the system introduced via the escort probabilities [TU] 
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may be recast in the following form: 
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(2.9) 
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Parallel to (12. 5p the above integral may be completed, and this, in turn, produces an 
implicit equation for the internal energy: 



9 /^Gs) 
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Substituting ( 12. 6ft in ( 12.101) the internal energy is found to be proportional to the trace c: 

Uf) = C N P' 1 . (2.11) 

Employing (12. 5ft . ( 12.71) and (12.111) we obtain an explicit expression for the generalized 
partition function: 



Zj 3) (/?) = [exp„(A0n /3 _iV£(1) , 



(fc)^(fc)- 



(2.12) 



The structure of $ given in ( 12.61) makes it evident that the generalized partition function 
( 12.12ft for A" rotators has singularities at the values of q given by q = (n + l)/n for 
n = 1,2,...,N. The number of singularities equals the number of degrees of freedom of 
the system. As the number of rotators A" increases the singularities tend to accumulate 
at q = 1, the limiting value where the statistical mechanics becomes extensive. The 
generalized partition function is, therefore, well-defined in the interval < q < 1 + 1/N, 
where the upper limit is characterized by the total number of the rotators. The above 
discussion makes it apparent that in the context of nonextensive thermodynamics the size 
A" of the system is of key importance in selecting the admissible range of the nonextensivity 
parameter q. In contrast to our result the authors of Ref. [TJ] considered the generalized 
partition function at the high temperature limit for the parameter range 1 < q < 2. 
Employing the results (12.71) . ( 12. lip and ( 12.121) the internal energy 
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and the specific heat of the system may be computed: 
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The specific heat of the rigid rotators in the high temperature limit, in contrast to its 
classical value, becomes a temperature dependent quantity for nonextensive statistics. In 
the extensive q — > 1 limit we recover the classical Boltzmann-Gibbs statistics for these 
thermodynamics quantities: U q —>■ N ksT, C q —>■ N ks- 



To obtain the internal energy (12.131) and the specific heat (12.141) we may also employ the 
method discussed in [10] , where the physical quantities evaluated in the second constraint 
approach are transformed to those in the third constraint approach by introducing a 
fictitious temperature f3' that is useful for calculation. The ensemble probabilities in the 
second constraint and the third constraint approaches are now linked as 



(2), 
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In other words, the second constraint approach is employed as a convenient tool at the 
intermediate stages of calculation, and finally results are translated to the third constraint 
approach via certain transformations. The general transformation rule for the tempera- 
ture reads 



l - 
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In the present example of the rigid rotators at the high temperature limit this relation 
assumes the form 
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The specific heat and the internal energy evaluated using this technique is identical to 
the results ( 12. 131) and (12. 14ft where the third constraint method has been used ab initio. 



Ill Rigid rotators in the low temperature limit 

In the example of N rigid rotators at the low temperature limit we obtain the thermody- 
namic quantities perturbatively by treating (1 — q) as the series parameter. This process 
involves use of the second constraint approach as an intermediate step. The results cor- 
responding to the third constraint approach are then extracted [10J via transformations 
(I2.15P and (I2.16P linking the ensemble probabilities in these two approaches. The partition 
function of iV rigid rotator in the second constraint setting reads 

oo 

Zf\0) = Trexp,(-/3 C) = J> [l-(l-q)$C]^, exp q (X) = (l+(l-g) X)^ (3.1) 

and the pertinent ensemble probability is given by 

p«C9) = [1 - (1 - q)PC)^ (Zf\p))-\ (3.2) 

In the low temperature regime the k = and U = 1 quantum levels generate the dominant 
contribution to the partition function. 
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Using an infinite product expansion of the g-exponential 

(1 -qf- 1 



exp q (-fi£) = exp ( - g !LJ£lpc k> ) 



(3.3) 



we express the partition function (13. ip as a pertubative series in (1 — q), and retain terms 
upto second order: 



h=o,i 

Recasting the above series as 



l-(l-q)tc 2 + (l-q) 2 (^£ 4 -^£ 



exp(-/3£). (3.4) 



B 2 d 2 



(3 s d 3 (3 A <9 4 
+ 



3 dp 3 8 df3 4 



]T Pexp(-/3£) (3.5) 
h=0,l 



2 dp 2 

we obtain the series expansion 

Z?\P) = [f(P)] N [l - (1 - q)ft + (1 - qf{ft + ft + ft + ft)], 
where the coefficients read 

ft = 6N(3 2 £{(3) + 18N(N - l)/3 2 [£(/?) ] 2 , 

ft = (6Np 4 -8Np 3 )£(f3), 

ft = (12QN(N- l)/3 4 - 72N(N- l)f3 3 )[£((3)} 2 , 

ft = (324iV(A^-l)(iV-2) ( 5 4 -72A^(A^-l)(iV-2)/5 3 )[^(/?)] 3 , 

ft = 1Q2N(N- l)(N-2)(N-3)$ i [S(j3)] i , 

£((3) = exp(-2/3)[/(/3)]- 1 , f(J3) = 1 + 3exp(-2/3). 

The defining property of the internal energy in the second constraint framework 

and the series (13. 6p allows us to extract the perturbative result 

U®(J3)=Gi + (l-q) (£ 2 + £ 3 + £ 4 ), 
where the coefficients may be listed as 
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The thermodynamic quantities in the third constraint scenario materialize on transforming 
(3 —>■ ft' a la (12.161) . The key quantity k enacting this transformation may be factorized as 

« = E^ (2) )" = M (^G 9 ))"'. XW) = £ 2>[1 - (1 - q)PC} q,1 - q - (3-10) 

i k=0,l 

In the present case we evaluate the sum x{0) perturbatively and retain terms upto order 
(1 — q) 2 . Recasting SIS db differential series 



d ft 2 d 2 
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we obtain x(P) as 

Xifi) = ZWtf) + [f(P)] N [(l - q)6N $£(p) + (1 - q) 2 12N{0 2 - p 3 )£(P) 

+3(N - 1) 2 - 3(3 3 )[£(P)] 2 - 9{N - l)(N - 2) p 3 [£{P)} 3 }}. (3.12) 

The relations (13.101) . (13.121) and the partition function (13. 6ft now yield « upto the order 
(1-q) 2 : 

k = 1 + (1 - q)Vi + (1 - g) 2 (P 2 + V 3 + Vt), (3.13) 
where the coefficients read 

P 1 = Nhyf(J3)+6N$£(l3), 

N 2 . . . 

V 2 = — [Inf(P)} 2 + (6NP 2 - 12NP 3 + 6N 2 p In /(/?)) £(/?), 

p 3 = (18iV(iV - l)/3 2 + (108iV - 72iV 2 )/3 3 )[£(/?)] 2 , 

P 4 = 216iV(iV-l)/3 3 [£(/3)] 3 . (3.14) 

The transformation (12.161) may now be produced as a perturbative series: 

P = P' + (l-q) P'g(p') + (1 - g) 2 ^(^), (3.15) 

where the coefficients read 

g(P') = N\nf(P') + 12NP'£(p'), 

h(p') = 432N(N -l)P' 2 [£(p')} 3 + 18NP >2 (5N -3 + (12-8N)/3')[£(p')] 2 

. 2 * 7V 2 

+6iV(3/3' -4P + NP') £(p') + — [\nf(p')} 2 . (3.16) 

To proceed with our computation, we, however, need to invert the series (I3.15p . This, 
upto the order (1 — q) 2 is given as 

= p[l -(1-q) g(P) + (1 - qf (Pg(P)g'(P) + [g(P)] 2 ~ KP))\- (3-17) 
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Insofar as the transformation (3 — > /3' a temperature dependent process we can not use 
the same construction of the transformation in the distinct regimes of high and low 
temperatures. This was, however, done in the earlier work on the rigid rotator [14J. On 
this we disagree with them. 

Employing the key properties (12. 7p and (12.151) we now extract the generalized partition 
function in the setting of the third constraint: 

\Zf (/3)p-«) = k + (1 - q) 2 l2N 2 p 2 [12p£(P) + 1b /(/?)] {£((3) - 3[£{P)] 2 }. (3.18) 

The defining thermodynamic relationships 

El - (319) 

C® = (3-20) 

now allow us to produce the specific heat at the low temperature limit in the context of 
the third constraint: 

Cf = C BG - (1 - q)[C x + C 2 + C 3 + C 4 ], (3.21) 
where the familiar classical Boltzmann-Gibbs extensive term Cbg reads 

C BG = 12Nkp 2 £{P) - 36Nkp 2 [£{p)] 2 . 

The nonextensive terms in (13. 21 p at the order of (1 — q) are listed below: 

& = 12Nk[N0 2 -2p 3 )lnf{p) + p 2 - 4/3 3 + 2p A ]£{p), 

C 2 = 36Nk[Nlnf{P){6p 3 -p 2 ) + {12 + 2N)p 3 -{U + 8N)p A ][£{P)} 2 , 

C 3 = 216Nk[p 3 + 2p 3 In f{p) + [ION - 12)/3 4 - 4/3 3 ] [£(/?)] 3 , 

C 4 = -3888N(N + l)k p 4 [£(p)]\ (3.22) 

In the classical q — > 1 limit the specific heat (I3.2ip of N rotators reduces, expectedly, to 
its Boltzmann-Gibbs value. In the N = 1 rotator case, we, unfortunately, fail to notice 
this occuring in |14j . 

At the low temperature limit the leading term in internal energy may now be computed 
via integrations in (13. 20H . We quote the result: 

Uf ] = U BG + (l-q) (Ut + U 2 +U 3 +U 4 + Ik), (3.23) 

where the coefficients read 

U BG = 6Ntt R exp(-2P) [1 - 3exp(-2/3) + 6 exp(-4/3)], 
^ = 12Nn R exp(-2P) 0-P% 

U 2 = Nn R exp(-i$) [1 + 36(N - 2)p + 18(8N + l)p 2 ], 

U 3 = 2Nn R exp(-Q(3) [(2 - N) + (11 - 20AQ6/3 - (12 + 21)18/3], 

U A = 81NnR eX p(-8/3) [(6JV + 3/2) + (UN + 3/2)$ + (36N + 42)/3 2 ], 



W 5 = 2916 ^" fl exp(-10/3) [(1 + 10P)(9N + 4) - 200(iV + l)p 2 }. (3.24) 
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IV Gas molecules with translational and rotational 
degrees of freedom 

Physically another aspect of nonextensivity may be of interest. Most often systems incor- 
porate multiple degrees of freedom. In the nonextensive scenario this may indicate that 
the thermodynamic quantities such as specific heat are likely to have temperature depen- 
dent terms of mixed origin. This terms are generated due to the coexistence of different 
degrees of freedom. Another way of visualizing this is that nonextensivity of statistics 
induces effective interaction terms between two separate component of the microscopic 
Hamiltonian of the system. In contrast to this, different degrees of freedom remain dis- 
joint in the usual extensive statistical mechanics. As a prototype of this idea we consider 
here an ideal diatomic gas endowed with both the translational and the rotational degree 
of freedom. For the purpose of simplification we assume the molecules to be of large mass 
so that we may disregard the vibrational modes. We also restrict our considerations here 
to the high temperature limit. 

The energy of a diatomic molecule is given by 



^ 2m 21 v ; ' 

i=l 



0,1, 



(4.1) 



The generalized partition function in the third constraint is given by 



1 



N\ h 3N 



J2 V 



1 - (1 - q) 



2 £ rr 

c \2m Ur 
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A? 



Hd 3 q4 3 Pi, (4.2) 
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where P 2 = J^iPi- Following Sec. [TTJ here also we, in the high temperature limit, replace 
the summation on the rotational levels by an integral: 



h=0 Jp 



V 



P 2 



c \2m 



U, 
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Y[dlid 3 Pi 



(4.3) 
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where the numerical factor reads V = V N /N\ h 3N . The partition function may now be 
exactly evaluated as 



-5N/2 



exp„ 







H5/2) 



G p~ N /3 



-3N/2 



(4.4) 



where G is a ratio of two gamma functions and (3 is a scaled temperature: 

-l 



q 



q 



l 



q 5N\ 



2-nm 



A similar evaluation of the internal energy leads us to 



(4.5) 



5N V c^ +1 



pz, 
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exp. 
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(4.6) 
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Fortunately the ratio of two implicit equations (14. 4p and (14.61) assumes a simple form 

5Nc 



£/(3) 

q 2(3 



(4.7) 



On substituting (14 .7p in (14. 4p and (14.61) we evaluate the generalized partition function and 
the internal energy explicitly: 



^ [exp g (5iV/2)f 5 / 2 /3 
5N 

w 



-N£, 



^(1-5) A J/J ^(g-l)AT£ (5/2) ^(q-l)3^£ {5/2) ^ 



(4.t 



The specific heat of the diatomic gas in the high temperature limit is directly produced 
via the internal energy in (14.81) 



CV>{/3) 



d 



Uf\(3) 



5Nk 



'(5/2) Aj?£ /2) ^ £(5/2) pl*- 1 )™**/* ^(?-i)iV£ C6/2) _ (49) 



dT ~ q 

The above specific heat of the diatomic gas can also be verified by employing the /3 — *■ (3' 
transformation procedure discussed and used in Sec. IIII1 The relevant transform in the 
present case is quoted below: 
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J' 



v /2 n£ g (i - g r 5iv / 2 



( 1 -9)-C(5/2) 



-(5/2) 



(4-10) 



Expanding the specific heat (14.91) of the gas molecules having both translational and 
rotational degrees of freedom in the parameter (1 — q), we retain terms upto the order of 
(1-g) 2 : 



25N 2 5N\ 

~8 rjj' 



(4.11) 



where = InV — (3N/2)\n (3 — Nln/3 + 5N. For the purpose of comparison we now 
consider an expansion of the translational specific heat [TT] in the parameter (1 — q) while 
retaining terms upto order (1 — q) 2 \ 

9N 2 3iVY 



c, 



3N 
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1 + (1 - q)Et 



q) 



CT2 

2 / 



3N_ 
2 + ~2~~* 



- 



4 J 
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where H$ = InV — (3N/2) ln/3 + 3N. Applying the same expansion in the rotational 
specific heat in the high temperature case (12.141) . we obtain 



C r = N k 



2 N 3 N 2 

l- q )~ r + (l- q f[=L + N~ r -- - 



(4.13) 



where S r = — Nln$ + 2N. Summarizing the above expansion in the region q —>■ 1 we 
observe that the specific heat of the gas molecules with both translational and rotational 
degree of freedom may be written as the following sum: 



C d = C t + C r + M t + M r + M t , r , 



(4.14) 
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where the additional terms on the rhs read 



M r 





3(ln/3) 2 -531n/3 , 



1 + 2AN In 
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2 k + (l-q) 2 N 2 k 
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Evidently the specific heat Cd of the diatomic molecules is not a linear sum of the trans- 
lational specific heat C t and the rotational specific heat C r . The additional terms owing 
their origin to nonextensivity may be divided into three classes. The terms M t and M r in 
(14. 141) contain the effects of purely translational and purely rotational degrees of freedom, 
respectively. Qualitatively, these terms may be described as a kind of nonlinear effect of 
mixing of degrees of freedom. On the other hand, the term Mt, r incorporate the effect of 
both the degrees of freedom. These cross terms indicate the existence of some sort of effec- 
tive interaction between individual degrees of freedom, such as the translational and the 
rotational examples discussed here. These effective interactions are temperature depen- 
dent and may be thought of as induced by the nonextensive statistics. This qualitatively 
distinct feature is absent in Boltzmann-Gibbs statistics. 



We have studied the specific heat of iV rigid rotators in the high temperature limit. 
Using the disentanglement of the q exponential in an infinite product series of ordinary 
exponentials, we have developed a pertubative method to calculate the specific heat upto 
any particular order in the nonextensivity parameter (1 — q). This procedure is applied to 
find the specific heat of N rotators in the low temperature limit. Finally we calculate the 
specific heat of a diatomic gas with a combination of the translational and the rotational 
degrees of freedom. The q — > 1 limiting value of the specific heat specific heat leads 
us to conclude that there is a mixing between the different kinds of degrees of freedom 
due to nonextensive statistics. This mixing of independent degrees of freedom may be 
visualized as an effective interaction between them owing its origin due to nonextensivity 
of statistics. A study of these problems in the context of two parametric entropies [ToTfTT] 
should be worth pursuing. 
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V Concluding Remarks 
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